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Abstract. The reductions of the multi-component nonlinear Schrodinger 
(MNLS) type models related to C.I and D.III type symmetric spaces are stud- 
ied. We pay special attention to the MNLS related to the sp(4), so(10) and 
so(12) Lie algebras. The MNLS related to sp(4) is a three-component MNLS 
which finds applications to Bose-Einstein condensates. The MNLS related to 
so(12) and so(10) Lie algebras after convenient Z§ or Z4 reductions reduce 
to three and four-component MNLS showing new types of -interactions 
that are integrable. We briefly explain how these new types of MNLS can 
be integrated by the inverse scattering method. The spectral properties of 
the Lax operators L and the corresponding recursion operator A are outlined. 
Applications to spinor model of Bose-Einstein condensates (BEC's) are dis- 
cussed. 

1. Introduction 



When spinor BEC's are trapped in magnetic potential, the spin degree of freedom 
is frozen. However, in the condensate trapped by an optical potential, the spin 
is free. We consider BEC's of alcali atoms in the F = 1 hyperfine state, elon- 
gated in x direction and confined in the transverse directions y, z by purely optical 
means. Then, in the absence of external magnetic fields is characterized by the 
magnetic quantum number which has three allowed values mp = 1,0, — l.Thus 
the assembly of atoms in the F = 1 hyperfine state can be described by a normal- 
ized spinor wave function <J?(x, t) = ($i(x, t), $o(x, t), <3?__i(:r, t)) T whose com- 
ponents are labelled by the values of tuf- In short the dynamics of such BEC's is 
described by a three-component Gross-Pitaevskii (GP) system of equations. In the 
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one-dimensional approximation described above the GP system goes into the fol- 
lowing three-component nonlinear Schrodinger (MNLS) equation in (ID) x-space 
[1]: 

id t $i + dl®! + 2(|$!| 2 + 2|$ | 2 )$i + 2$L 1 $g = 0, 
id t $o + d 2 x $ Q + 2(|$_i| 2 + |$ | 2 + |$i| 2 )$o + 2$5$i$_i = 0, (1) 
+ a 2 $_i + 2(|$„!| 2 + 2|$ | 2 )$-i + 2$I$g = 0. 

This model is integrable by means of inverse scattering transform method [1]. It 
also allows an exact description of the dynamics and interaction of bright solitons 
with spin degrees of freedom. Matter-wave solitons are expected to be useful in 
atom laser, atom interferometry and coherent atom transport. It could contribute 
to the realization of quantum information processing or computation, as a part of 
new field of atom optics. Lax pairs and geometric interpretation of the model 
(1) are given in [2]. Darboux transformation for this special integrable model is 
developed in [3]. We will show that the system (1) is related to symmetric space 
C.I ~ Sp(4)/U(2) with canonical ^-reduction and has natural Lie algebraic 
interpretation. 

The applications of the differential geometric and Lie algebraic methods to soli- 
ton type equations lead to the discovery of close relationship between the MNLS 
equations and the symmetric spaces [2]. It was shown that these MNLS systems 
have Lax representation with the generalized Zakharov-Shabat system as the Lax 
operator: 

L*l>(x, t, A) = i^p + (Q(x, t) - \J)iP(x, t, A) = 0, (2) 

where J is a constant element of the Cartan subalgebra h C of the simple Lie 
algebra g and Q(x,t) = [J,Q(x,t)]. In other words Q(x,t) belongs to the co- 
adjoint orbit Mj of g passing through J. 

The choice of J determines the dimension of M.j which can be viewed as the 
phase space of the relevant nonlinear evolution equations (NLEE). It is equal to 
the number of roots of g such that a(J) / 0. Taking into account that if a is a 
root, then —a is also a root of g; therefore dim M. j is always even. 

We concentrate on those most degenerate choices for J for which ad j has just two 
non-vanishing eigenvalues ±2a; in this case J 2 = a 2 11. Such choices of J are com- 
patible with several types of symmetric spaces: A. Ill ~ SU(p + q)/S(U(p) <g> 
U(q)), C.I ~ Sp(2p)/U(p) and D.III ~ SO(2p)/U(p) [2, 4]. The classification 
of the symmetric spaces related to a given simple Lie algebra g is directly related 
to the classification of the Cartan involutions (i.e. to the classifications of the real 
forms) that the algebra admits. For more details see e.g. [4, 5]. 
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The interpretation of the ISM as a generalized Fourier transforms and the expan- 
sion over the so-called "squared" solutions (see [6, 7] for regular and [8, 9, 10] 
for non-regular J) allow one to study all the fundamental properties of the corre- 
sponding NLEE's. These include: i) the description of the class of NLEE related 
to a given Lax operator L(A) and solvable by the ISM; ii) derivation of the infinite 
family of integrals of motion; and iii) their hierarchy of Hamiltonian structures. 

The degeneracy of J means that the subalgebra qj C of elements commuting 
with J (i.e., the kernel of the operator ad j) is non-commutative which makes more 
difficult the derivation of the fundamental analytic solutions (FAS) of the Lax oper- 
ator (2) and the construction of the corresponding (generating) recursion operator 
A. Here we continue our studies in [7, 8] finding new algebraic reductions of 
MNLS equations related to C.I and D.III type symmetric spaces. Some of them 
like eq. (1) find applications to Bose-Einstein condensates and nonlinear optics. 
The derived reduced MNLS system seem to be new to the best of our knowledge. 

The present article is organized as follows: In section 2 we give some preliminaries 
about the simple Lie algebras and the general form of the MNLS models and the 
relevant recursion operators. Section 3 is devoted to the spectral properties of the 
Lax operator L. In section 4 we discuss the Hamiltonian properties of the MNLS 
systems using the classical i?-matrix method. In section 5 we apply the approach 
in [11, 9] and derive the dressing factors and the soliton solutions for symmetric 
spaces related to D5 and Y)q Lie algebras. Section 6 is devoted to the analysis 
of the reductions of the MNLS equations by applying the reduction group [12] 
method. The relation between reductions and scattering data for L operator are 
outlined in section 7. 



2. Preliminaries. 

2.1. Simple Lie algebras. 

Here we fix the notations and the normalization conditions for the Cartan-Weyl 
generators {h k , E a } of q (r = rankg) with root system A. We introduce h k 6 h, 
k = 1, . . . , r as the Cartan elements dual to the orthonormal basis {e^} in the root 
space E r and the Weyl generators E a , a G A. Their commutation relations 




r 



[h k ,E, 



'a 



(a,e k )E ( 



k=i 




for a + G A 
fora + /3 Au{0}. 



(3) 
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Here a = J2k=i a k^k is a r-dimensional vector dual to J G fj and (-, •) is the scalar 
product in E r . The normalization of the basis is determined by: 

E. a = E%, {E_ a ,E a ) = -^—, N- a> _p = -N a! p,, (4) 

where N a $ = ±(p + 1) and the integer p > is such that a + s/3 G A for all 
s = 1, . . . ,p, a+ (p + l)(3 A and (-, •) is the Killing form of g, see [4]. The root 
system A of g is invariant with respect to the group W 5 of Weyl reflections S a , 

b a y = y- -, T-a, a G A. (5) 

(a, a) 

With each reflection S a one can relate an internal automorphism of the algebra 
Ad A a £ Auto which act in a natural way on the Cartan-Weyl basis, namely: 

S a (Hp) = AuHpAa 1 = Hp>, (3' = S a p, 

S a {Ep) = AaEpA' 1 = n aj pEpi, n a> p = ±1. (6) 

Since 5„ = 11 we must have A\ = ±1. 

As we already mentioned in the Introduction the MNLS equations correspond to 
Lax operator (2) with non-regular (constant) Cartan elements J € f). If J is a regu- 
lar element of the Cartan subalgebra of g then ad j has as many different eigenval- 
ues as is the number of the roots of the algebra and they are given by aj = (Xj(J), 
a.j G A. Such J's can be used to introduce ordering in the root system by assum- 
ing that a > if a (J) > 0. In what follows we will assume that all roots for which 
a(J) > are positive. 

Obviously we can consider the eigensubspaces of ad j as grading of the algebra g. 
In what follows we will consider symmetric spaces related to maximally degener- 
ated J, i.e. ad j has only two non-vanishing eigenvalues ±2a. Then g is split into 
a direct sum of the subalgebra go and the linear subspaces g±: 

S = So g+ fl_, g± = I.e. {X ±j I [J, X±] = ±2aX±} 

The subalgebra go contains the Cartan subalgebra t) and also all root vectors E± a 6 
g corresponding to the roots a such that (a, a) = 0. The root system A is split into 
subsets of roots A = #o U 9 + U (— 9 + ), where: 

9 = {a g A | a( J) = 0} , 0+ = {a G A | a(J) = a > 0} . (7) 

We can use the gauge transformation commuting with J to simplify Q; in particular 
we can remove all components of Q in go; effectively this means that our Q(x,t) = 
Q+(x, t)+Q-(x, t) G g+Ug_ can be viewed as a local coordinate in the co-adjoint 
orbit Mj ~ g\g : 



Q+(x,t)= q a {x,t)E a , Q_(x,t)= ^2 p a (x,t)E a . (8) 
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Obviously Q± £ q± and: 

adjQ= [J,Q] =2a(Q+-Q_), (ad j^Q = ±-{Q+ - Q_); (9) 

la 

besides [E a , Ep] = for any pair of roots a, j3 G 0+. This simplifies solving the 
recursion relations and the explicit calculation of the recursion operator A. 



2.2. Lax representation of the MNLS type models 

The operator (2) together with the corresponding operator M(A): 

M(X)ip = (i^- t - [Q, ad fQ] + 2iad7 1 Q a; + 2XQ - 2X 2 jj ^{x, t, A) = 0, 

(10) 

where Q = Q(x,t), provide the Lax representation for the MNLS type systems. 
The compatibility condition [L(A),M(A)] = of (2) and (10) gives the general 
form of the MNLS equations on symmetric spaces: 



J ' dt 



d 2 Q 



+ ^}-2a 2 [ad7 1 Q,[ad7 1 Q,Q]] = 0. 



(11) 



Following [13] one can consider more general M-operators of the form: 

HvJ> ( N \ 

M(A)tt = i— + ]T v "( x > ^ *> A ) = 0, /(A) = Jim^ V(x, t, A). 

\fe=i / 

(12) 

The Lax representation [L(A), M(A)] = leads to a recurrent relations between 

V k (x,t) = V k i + V k d 

V k t +1 (x,t) = TTj(V k+1 ) = A±V k i (x,t)- a dj 1 [C k ,Q(x,t)], (13) 
V?(x,t) = (i-irj)(V k ) = C k + i f dy[Q(y,t),V k { (y,t)}, k = l,...,N; 

J ±oo 

where ttj = ad J 1 o ad j and C k = (11 — ir,j)C k are block-diagonal integration 
constants, for details see, e.g. [13, 2]. These relations are resolved by the recursion 
operators (14): 

adj f.dZ 



K+Z 



4a 2 V 1 dx +1 



Q(x), f dy[Q(y),Z(y)} 

J±oo 



(14) 



where we assume that Z = ttjZ £ M j. As a result we obtain that the class of 
(generically nonlocal) NLEE solvable by the ISM have the form: 



C fe ,ad7 1 Q(x,t) =0, f(X) 



/+(A) 
/-(A) 



(15) 
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where /(A) = Y,k=o Ck^ N k determines their dispersion law. The NLEE (15) 
become local if /(A) = /o(A) J, where /o(A) is a scalar function. In particular, if 
/(A) = -2A 2 J we get the MNLS eq. (1 1). 



2.3. Basic physical example: C.I type symmetric space Sp(2p)/U(p). 

We choose Q = C 2 = sp(4) algebra; it has 2 simple roots, namely a\ = e\ — e 2 , 
a>2 = 2e 2 . We fix up the Cartan element 

J = diag (a, a, — a, — a), J 2 = a 2 ll. (16) 

Then the corresponding potential Q(x,t) (8) takes the form: 



Q(x,t) 



/ 9 12 ?i ^ 

q 2 -qi2 

P12 P2 

V pi -P12 J 



(17) 



Imposing the involution p& = q^, k = 1, 2 and P12 = g* 2 we obtain the following 3 
component MNLS system for the independent fields (712(2;, t), (71 (x, i) and q 2 (x, t): 

. dqu d 2 qi2 ,, 12 , 1 12 , 1 i2\ * n 

iQ^r— + „ n + 2gi 2 (|gi 2 | + \qi\ + ml ) ~ 2q 1 q 2 q 12 = 0, 



dt ' dx 2 

^ + ^ + 29l(kl ' 2 + 2|ft2|2) " 2 ^ 2 = ° 
dt dx 2 



(18) 



+ 2 (/2 (| (/2 | 2 + 2|( Zl 2| 2 )-2( Z 2 2(h * = 0, 



If we identify the physical quantities of the system (1) with 

$0 = 312, $1 = -qi, $-1 = Q2- (19) 

then eq. (18) will coincide with (1). The system can be written in a Hamiltonian 
form by introducing the Poisson brackets: 

{qj(x),Pk(y)} = 2iS kj S(x - y), {q 12 (x),p 12 (y)} = iS(x-y), (20) 
and the Hamiltonian: 



H 

#kin 



int ! 



dx 



\ dx dx 2 V dx dx dx dx 



(21) 



= ~ /" dx ((|d>o| 2 + |*i| 2 ) 2 + (|$o| 2 + |<&-i| 2 ) 2 ) 

- - / dx (|$o$li + ^i^| 2 ) • 
The soliton solutions of the sp(4) MNLS (1) were derived independently in [1, 14]. 
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2.4. D. Ill type symmetric space 

We choose g = Dq ~ so(12); it has 6 simple roots, namely a\ = e\ — e 2 , 
a 2 = e 2 - e 3 , a 3 = e 3 - e 4 , q 4 = e 4 - e 5 , a 5 = e 5 - e 6 and a 6 = e 5 + e 6 . We 
fix up the Cartan element 

J = diag(a,a,a,a,a,a,—a,—a,—a,—a,—a,—a), J 2 = a 2 ll (22) 

which means that the subset 6 + = {e^ + ej} with 1 < i < j < 6. Then the 
corresponding potential Q(x, t) (8) takes the form: 
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(23) 



Here by qij (x, t) and pij (x, t) where i, j belong to the set of indices J = { (ij) | 1 < 
i < j < 6} we denote the coefficients of the generators E a and E_ a with 
a = ei + ej. Then the generic NLEE (11) becomes a system of 30 equations. 
Imposing the natural involution Q = Q\ i.e. = q*j we obtain a MNLS for the 
15 independent functions qij(x, t). 

Similarly we consider also the D5 ~ so(10) algebra. The corresponding Q(x,t) 
will be a 10 x 10 matrix-valued function which can be obtained from the equation 
for (23) by removing the first and the last rows and columns. The generic NLEE 
will be a system of 20 equations, the involution Q = will reduce them to a 
10-component MNLS. 

3. Spectral data and generalized exponents 

Here we will start with a brief sketch of the direct scattering problem for (2). It is 
based on the Jost solutions [15, 16] defined by their asymptotics 

lim ?/;(x,t,A)e iAJa; = H, lim <f>(x, t, X)e iXJx = 11, 
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and the scattering matrix T(t, A) = (Y>(x, t, A)) 1 (f){x, t, A) and its inverse T(A, t): 

T(f \\ _ ( A) -b~{t, A) \ f _ f c-(t, A) d~(t, A) \ 
T ^' A )- [b+(t,X) a-(t,\) )' { ' >~{-d + (t,\) c+(t,\) ) A2A) 

where a ± (t, A) and b ± (t, A) are r x r block matrices. The blocks a ± , b ± , and 
satisfy a number of relations coming from the fact that T(A)T(A) = 1; for 
example 

a + (A)c _ (A) + b-{\)d + {\) = 11, a+(A)cT(A) - 6~(A)c + (A) = 0, (25) 
etc. 

The fundamental analytic solutions (FAS) x^x,*, A) of L(A) are analytic func- 
tions of A for Im A ^ and are related to the Jost solutions by: 

X ± (x, t, A) = 4>(x, t, X)Sf(t, A) = 1>(x, t, X)Tj(t, A). (26) 

Here Sf, Tf upper- and lower- block- triangular matrices: 



a-(t,X) J ' J ^ ' > ~ \ b+(t,X) 11 

satisfy lf(t, X)Sf(t, A) = T(t, A) and can be viewed as the factors of a gener- 
alized Gauss decompositions of T(t, A) [8]. If Q(x,t) evolves according to (15) 
then 

i_ dT + ^ (A)6±(t ' A) " 6±(t ' A)/±(A) = °' i_ dT + [/±(A) ' a±(t ' A)] = °- 

(29) 

On the real axis in the complex A-plane both FAS x^^-, t, A) are linearly depen- 
dent: 

X + (x,t,X) = x -(x,t,X)G (t,X) (30) 

and Go(t, A) can be considered as a minimal set of scattering data in the case of 
absence of discrete eigenvalues for the Lax operator (2), see [10]. 
The mapping between the potential of the Lax operator and the scattering data is 
based on the Wronskian relations [17]. As an example of one we write down: 

i f 00 

bt(t,\) = -z / dx([Q(x,t), J]x + (x,t,X)E a ( X + (x,t,X))- 1 ), 

£ J -co 

where E a is the root vector corresponding to the root a £ 0+. Thus the "squared" 
solutions that appeared first in [13, 6] were later generalized in [7, 8] to Lax oper- 
ators of the type (2) as follows: 

et(x, t, A) = it j ( X ± (x, t, X)E a ( x ± (x, t, A))" 1 ) . (31) 
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They can also be viewed as natural generalizations of the usual exponentials and 
their completeness relations in Ai j [7, 8] provide us the spectral decompositions 
of the recursion operators A± for which e± a {x, A) are eigenf unctions: 

A+e± a (x, A) = \e± a (x, A), A_e| Q (x, A) = Ae| a (x, A), a G 6 +.(32) 

The (generating) recursion operators A+ appeared first in the AKNS -approach [13] 
as a tool to generate the class of all M-operators as well as the NLEE related to the 
given Lax operator. Next I. M. Gel'fand and L. A. Dickey [18] discovered that the 
class of these M-operators is contained in the diagonal of the resolvent of L. The 
kernel of the resolvent of L can be explicitly defined in terms of the fundamental 
analytic solutions x ± ( x i A) of (2), see [19, 7, 8]. 

4. Hamiltonian properties of the MNLS models 

It is well known that the MNLS equations possess hierarchies of Hamiltonian struc- 
tures. The phase space M. j of the MNLS equations is the co-adjoint orbit of the 
~ D r determined by J; in addition we assume that the matrix elements of 
Q(x, t) are smooth functions tending to zero fast enough for \x\ — ► oo. 
On the D. Ill-type symmetric spaces the Hamiltonian of (11) is given by: 

#m nls = « dx |2(ad7 1 Q a; ,ad7 1 Q a; ) + i([ad7 1 Q,g],[ad7 1 Q,Q])|. 

(33) 

Direct calculation shows that this Hamiltonian is proportional to the third coef- 
ficient I3 of the expansion of the generating functional of the principal series of 
motion In det a + with respect to A: 

oo 

In det a + (t,X) = j^ifcA - *, 
k=l 

i.e. -f^MNLS = <^3- The equation of motion (11) that Q satisfies is generated by 
the canonical symplectic structure: 

^mnls = ~ 2ia J°° dx(sQ(x,t) AadJ^QCM)^ . 

The hierarchies of symplectic structures defined on Ai j are generated by the cor- 
responding recursion operators A± (14) and are given by the following families of 
compatible two-forms: 

^NLS = - 2ia j°° dx^5Q(x,t) AA k adj 1 5Q(x,t)^ , A=^(A + + A_). 

(k) 

The corresponding Hamiltonians for the higher MNLS are ff^ LS = 8i/3 + fe. 
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For /(A) = — 2A 2 J eq. (29) gives da ± /dt = and a ± (A) can be viewed as gen- 
erating functional of integrals of motion whose number r 2 is larger than the rank 
r of g. This is obviously due to the degeneracy of the dispersion law. For generic 
/(A) from (29) there follows that only functions of the eigenvalues of a ± (\) will 
be conserved. Indeed, it follows from the classical i?-matrix approach [15]. One 
of the definitions of the classical i?-matrix is based on the Lax representation of 
the MNLS (2). In particular, if U(x, A) has the form 

U(x,X) = Q(x,t)-\J 

and the matrix elements of Q(x, A) satisfy the Poisson brackets {•, •}: 

{Pa 2 (x, t), q ai (y, t)} = i 5 aia2 5{x - y) (34) 

then the classical i?-matrix can be defined through the relation 

{U(x, A) ® U(y, n)} = [R(x - y), U(x, A) ® II + II ® U(y, fx)}S(x - y) (35) 

Our system of equations (35) allows i?-matrix given by [2]: 

*(A-M) = ^P, 06) 

where 

?=(±h k ®h k+ ^f^A. (37) 

\fc=l aeA^^-aU 

Here h k are introduced in (3) and are properly normalized to (hi, h k ) = &ik- P 
is the second Casimir endomorphism of the algebra and posses special properties 
concerning its action onto the matrices of the corresponding group: 

P(A ®B) = {B® A)P, 

Using these properties of the P-matrix and the commutation relations of the alge- 
bra (3) we obtain 

[P,Q(x) II + II ®Q(x)] = 

[P, A J ® 11 + /xl ® J] = 2(A - fi){U(x, A) (g> U(x, n)} 

These relations are true for any Q(x) taking value in the algebra and the relation 
(35) seems most natural and its r.h.side does not contain Q(x, t). 
Let us now show, that the classical i?-matrix is a very effective tool for calculating 
the Poisson brackets between the matrix elements of T(A). It will be more con- 
venient here to consider periodic boundary conditions on the interval [— L, L], i.e. 
Q(x — L) = Q(x + L) and to introduce the fundamental solution T(x, y, A) [15]: 

dT(x,y,X) 



dx 



+ U (x, A)T(x, y, A) = 0, T(x,x,X) = 1L. (38) 
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Skipping the details we just formulate the following relation for the Poisson brack- 
ets between the matrix elements of T(x, y, A): 

T(x, y, A) ® T(x, y, = [R(X - n),T(x, y, A) T(x, y, /x)] (39) 

The corresponding monodromy matrix Tl(\) describes the transition from — L to 
L and 1l(A) = T(—L, L, A). The Poisson brackets between the matrix elements 
of Tl{\) follow directly from eq. (39) and are given by: 

T L {\) ® r L (/x)| = [R{\ - M ), T L (A) ® T L (/i)] . (40) 

We could also write the Poisson brackets between the matrix elements of the in- 
verse of the monodromy matrix Tl(A): 

{t l (A)®T l ( m )} = [Ti,(A) <8) Tl(h),R(\ — //)] . (41) 

An elementary consequence of this result is the involutivity of the integrals of 
motion 1^^ and from the principal series which appear in the expansions of: 

oo oo 

lndeta+(A) = ^/ L , fc A- fc , - lndet c£(A) = ]T J^A"*, (42) 

fc=l k=l 



oo oo 



lndetc+(A) = ^J L , fc A- fe , -lndeta^A) = £ J^A"*, (43) 

k=l k=l 

As a result of the reduction conditions (75) that we impose, the generating func- 
tional of the principal series of integrals of motion is only one, i.e. In det a£ (A). 
Another important property of the integrals 1^ and is their locality, i.e. their 
densities depend only on Q and its x-derivatives. 

The simplest consequence of the relations (40) and (41) is the involutivity of 1^ 
and Jl^- Indeed, taking the trace of both sides of (40) and (41) shows that 
{trT L (A),trT L (/i)} = and {trf L (A),trf L (^)} = 0. We can also multiply 
both sides of (40) and (41) by C ® C and then take the trace. This proves: 

{trT L (A)C,trT L (/x)C} = 0, {trf L (A)C , ,trf L (/i)c} = 0. (44) 

In particular, for C = II + J and C = 11 — J we get the involutivity of: 

{tra±(A),tra±(/i)| = 0, {tra£ (A),tr a£(p)} = 0, 
L J 1 J (45) 

{trc±(A),trc±(Ai)} = 0, {trc+(A),trc^(/i)} = 0, 

Equations (40) and (41) were derived for the typical representation of the 
corresponding group G, but they hold true also for any other finite-dimensional 
representation of G. Let us denote by the r-th fundamental representation of 
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G; then the element Tl(A) will have representation in V^'\ see [4]. In particular, if 
we consider equations (40) and (41) in the representation and sandwich them 
between the highest and lowest weight vectors in we get: 



{deta£(A),deta£(,u)} = 0, {det c£(X), det c£(/x)} = 0. 



(46) 



Since eq. (46) hold true for all values of A and /x we can insert into them the 
expansions (42) with the result: 



{lL,ki Il, p } = 0, {JL,k, Jl, p } = 0, 



k,p = 1,2, . . . . 



(47) 



Lets go back to our basic physical example from Section 2.3 associated with the 
algebra g = sp(4) and consider the monodromy matrix T^(A) in the typical repre- 
sentation. The matrix elements of the diagonal blocks of Tl a t,(\) for 1 < a < 2 
and 1 < b < 2 are denoted by Tl a b(ty = a ab accor ding to (24). Now Poisson 
brackets between the scattering data aJ(A) and a£(/t) are obtained in a straight- 
forwar way from: (40): 

a Lll( A Wl 2 (^)} 
a Lll( A )' a L2l(^)} 
a Lll( A )' a L22(^)} 
a L22( A )' a L12^)} 
a L22( A )> a L2l(^)} 
a L2l( A )' a L12(^)} 



(a+ n (A)a+ 12 (/x) - a+ 12 (A)a+ n (/x) 
1 — - (a+ 21 (A)a+ n (/x) - a+ n (A)a+ 21 (/x) 
(a+ 21 (A)a+ 12 (/x) - a+ 12 (A)a+ 21 (^) 





2(A- 




2(A- 




2(A- 




2(A- 




2(A- 
i 


m) 


2(A- 





(a£ 22 (A)a£ 21 (//) - a+ 21 (A)a+ 22 (^) 
(a+ n (A)a+ 22 (/x) - a+ 22 (A)a+ n (/x) 



a Lll( A )' a Lll(^)} = °> { a ^22( A )> a L22(^)} 
«^l 2 ( A )' a L12(^)} = 0. { a ^2l( A )' a L2l(^)} 





0. 



(48) 



Now that we know the Poisson brackets between the matrix elements of a + it is 
not difficult to extend the above result over its invariants: 

{tra+(A),a+^)} 



2 ;^3^[°£( a )>°l (/*)]«> 



for 



l,2j = l,2 



{deta+(A),a+. J .( M )} = 0, for i = 1,2 j = 1,2 



This somewhat more concrete analysis allows one to see that only functions of the 
eigenvalues of a£(A) produce integrals of motion in involution. 
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We are able to transfer these results also for the case of potentials belonging to the 
co-adjoint orbit of the g ~ D r determined by J and with zero boundary conditions 

{T(A, t) ® T(ji, t)} = [R(X - /i), T(A, t) ® T(n, t)}. (49) 

Here T(A, t) is the scattering matrix, obtained after taking the limit L — > oo in the 
corresponding monodromy matrix ?l(A). Indeed, let us multiply (39) by E(y, A)® 
E(y, fi) on the right and by E^ 1 (x, A) ® E" 1 ^, //) on the left, where E(x, A) = 
exp(— iA Jx) and take the limit for x — > oo, y — > — oo. Since: 

„io;(A— //) 

lim — = ±i7rJ(A - n), (50) 

x^ioo A — jJL 

we get: 

T(A) ® T(/x) J = i?+(A - /i)T(A) T(/x) - T(A) T(/x)i?_(A - /i), 

! / r x (51) 

#±(A - m) = 2 ^ A _ ^ ^Xj h k®h k + n 0J J ± i7T(5(A - n)Uu, 

where IToj and IIu are defined as follows 

n 0J = ( E <* ® + ® ^a) , 

^ (52) 

nu = ^ (E a ® - £_ a ® # Q ) . 

Analogously we prove that: 

i) there are integrals I k = Hindoo II^ and J p = lirn^oo Jl :P that are in involu- 
tion, i.e.: 

{I k ,I±} = {J k ,J p } = 0, 

for some positive values of k and p; ii) the principal series of integrals of motion 
Ik, generated by lndet a ± (t, A), i.e. the eigenvalues of a ± (A) and c^(X) produce 
integrals of motion in involution. So for MNLS we have extra integrals of mo- 
tion that are not all in involution. Indeed, lets denote the matrix elements of the 
scattering matrix T(A, t) according to (24) and multiply (51) by E ab ® E cd on the 
right, where (E ab )ij = 5 a i5 b j and take the trace of the elements in the first and 
in the second position of the tensor product. Thus we obtain the Poisson brackets 
between the scattering data: 

{T ba (X),T dc (p)} = tr (T(A) ® T(p) (E ab ® E dc R + - R^E ab ® E dc )) (53) 

Lets list the relations contained in the above equality and concerning the block- 
diagonal portions of the scattering matrix a + (X) and a - (A) 

{a±(A),a±>)} = (a± (A)a±» - a±(A)a± (/,)) (54) 
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{tr(a ± (A)) fe ,a±( A1 )} = l-L- [(a ± (A)) fe , o ± ( M ) 



The above result allows us to compute the Poisson brackets between the invariants 
of a ± (A) and their matrix elements: 

k_ 

2 A - n L v v " ' y ' 'lab (55) 
{trlna ± (A),a±(^)} = 

This analysis allows reveals that only the eigenvalues of a ± (A) and c =t (A) produce 
integrals of motion in involution. 

From (49) it follows that the first integrals Ik generated by lndet a ± (t,X) are 
in involution. Due to the special degenerate choice of the dispersion law /(A) = 
— 2A 2 J, any matrix elements of the blocks a ± (A) will generate integrals of motion, 
which however will not be in involution, see (54,55). The Hamiltonian for the 
MNLS models is proportional to 1% ', i.e. belongs to the principal series. If we 
choose a generic (i.e., non-degenerate) dispersion law then the Hamiltonian of the 
corresponding NLEE will not be in involution with I^. Such are the dispersion 
laws for the NLEE's that allow "boomeron" and "trappon" type solutions [20, 21, 
22]. This is the reason why their velocities become time-dependent. 

5. Dressing factors and soliton solutions 

n n Frt~\ m a PA Q x 



The main idea of the dressing method is starting from a FAS x%)( x i °f ^(^) 



with potential to construct a new singular solution ^ (x, A) of the Riemann- 
Hilbert Problem (30) with singularities located at prescribed positions A^. Then 
the new solutions xfi)( x i X) Wli l correspond to a potential of L(A) with two 
discrete eigenvalues Xf. It is related to the regular one by the dressing factors 

u(x, A): 

xfi)(x,X) = u(x, X)xf 0) {x,X)uZ 1 {X), it_ (A) = lim u(x,A). (56) 
If q ~ B r , D r the dressing factors take the form [23] : 

u(x,X) = 11 + (d(A) - l)Pi(x) + (c^(A) - l)P_i(x), P_i = Stff ST 1 , 

(57) 

where the rank 1 projector Pi(x) and the function ci(A) are given by: 

|n(x))(m(x)| A - A+ 

= <m(x)|n(x)> ' Cl(A) = A^F' 

K^)) = Xo{x,Xf)\n ), (m(x)\ = (m \x (x, Af). (58) 
| no) and (mo| are constant vectors and 

S = £(-l) fc+1 (P fc£ + ~k = 2r + 1 - k. 

k=i 
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Here E k n is an 2r x 2r matrix whose matrix elements are {E kn )ij — Si k 5 n j. Then 
the "dressed" potential have the form 

Q(i)OM) = Q(o){x,t)-{Xi-Xi)[J,p{x,t)}, p(x,t) = P 1 (x,t)-P- 1 (x,t). 

(59) 

where 

p(x,t) = -4- \J2hk(x,t)H ek + J2 (Pa(x,t)E a +P_ a {x,t)E_ a ) \ , 

\ m \ n ) y k=1 QgA+ j 

(m\n) = ^(nofcmofce^^+^^^' + no^mo^e- 2 ^^- 16 ^^'), 
k=l 

h k (x,t) = n 0k m 0k e 2a ^ x+Wa ^ t - n 0k m 0k e- 2a ^ x - 16a ^ t , 
and 

P a (x,t) = (no fe m 0s --(-l) s+/c n 0s m o5 )e- 2i ^^- 8ia ^-^) t , 
P- a (x,t) = (n 0s -m 0fc -(-l) s %m 0s )e 2i * 3;+8i ^-''?) t ; 
for a = e k + e s , and 

P a (x,t) = n 0k m 0s e 2a,/lX+Wa ^ Ult - {-l) s+k n^m 0k Q- 2a ^ x - ma ^ t , 
P- a (x,t) = n 0k m 0s& - 2a ^ x - ma ^ t - (-l) s+k n 0s m ke 2aulX+WafllUlt (60) 
for a = e k — e s , k < s. 

Let us consider now the purely solitonic case, i.e. Q^(x) = and 

xf )(x,t,\t) = exp(-iA^ Jx - 4iXi 2 Jt). 
Thus the 1 -soli ton solution is 

2ie^i (n 0j m ok - (-iy +k n ok m j) g-^ix-sia^-^)* 

Qjk = ^— — — 7 — — t , (61) 

^/ipi(p2ch \ 2av\x + + ^ ln^ J 

where (ij) € J and 

r r 

<Pi = ^2 n 0j m 0j , ip 2 = ^2 n oj m oji Xf = m± \v x . (62) 

3=1 3=1 

6. New reductions of MNLS equations 
6.1. The reduction group 

The reduction group Gr introduced by A. V. Mikhailov [12] provides a power- 
ful tool for constructing new integrable equations [11, 24, 23, 25, 26, 27] starting 
from known ones. It is a finite group which preserves the Lax representation, i.e. 
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it ensures that the reduction constrains are automatically compatible with the evo- 
lution. The main idea of the reduction group is to impose an invariance condition 
on the Lax operators (2) and (10). In particular this means that the dispersion law 
/mnls(A) = — 2A 2 J must also be compatible with the reduction group action. 
Here we consider two types of Gr reductions; like in [9] we will embed them as 
subgroup of W g : 

Type I: B^W {x,t, X*)B = U{x,t, X), B~ 1 JB = J, 

B- 1 V*{x,t,\*)B = V{x,t,\), (63) 
Type II: C^ 1 U*(x,t, X*)C = —U(x,t, A), C- 1 JC = -J J 

C- 1 y*(x,t,A*)C = -V(x,t,X), (64) 

where 

U(x,t, A) = Q(x,t) - A J 

and 

V(x,t,X) = -[Q, ad j l Q] + 2iadJ 1 Q ;c (x, t) + 2XQ(x, t) -2X 2 J. 
The automorphisms C and B must be of even order. 



6.2. Example of ^-reduction. 

Let us impose the following Z4 -reduction: 

B~\rf{x, t, X*))B = U(x, t, A), U(x, t, X) = Q(x, t) - A J 

B = W ei „ e2 ■ W e2 - e3 ■ W e3 - e4 , B 4 = 11, 

where w ei - ej are the Weyl reflection with respect to the roots — ej of the so(10)- 
algebra. Then B(J) = J and the corresponding reduced potential Q red G so (10) 
takes the form: 
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(65) 



NEW MULTI-COMPONENT NLS TYPE EQUATIONS ON SYMMETRIC SPACES 17 



Thus one derives the following 4 component MNLS system related to D5- algebra 
(here the independent fields are: qn(x, t), qu(x, t), qis(x, t) and q25(x, t): 

■ dqu d 2 qi2 |2 , ,2,1 ,2,1 i2\ 
m ~Qf + -^2~ + 2< lm\qi2\ +2|gi 4 | +\qi5\ + \Q25\ ) 

+2^i 4 (kl 5 | 2 + Ifel 2 ) + 2qf 4 q* 12 = 0, 

ia ^w + iS 1 + 2914(191412 + 2|912 ' 2 + |9is|2 + |(?2s|2) 

+2(?i 2 (ki 5 | 2 + k 25 | 2 ) + 2g 2 2 <zi4 = 0, (66) 

■ <9<?15 3 2 gi5 . . ,2 , ol 1 2 1 I |2 1 I |2\ 
la ~Qf + -^2~ + 2 <?15(2|gi5r + 2|925| + I912I + kl4|) 

+ 2^12915^14 + 2gi4^159l2 = 0, 

• ^25 . 9 2 q 2 5 , , i2 , ol 1 2 1 I i2 1 I |2\ 

+ -^3- + 2g 2 5(2|g 2 5| + 2|<?i 5 | + \q\2\ + \qu\ ) 

+ 2<?14<725<?*2 + 2gi2'?25<7l4 = °i 

The Hamiltonian of (66) is obtained from the general expression (33) by imposing 
the reduction constraints: 

H = i^kin + H int , Hki n = - [ dx + ' ' 



a 



00 



dx dx dx dx 



+ dqndqh + dmd^A 
dx dx dx dx J ' 

#int = / dx (|<?1 5 | 2 + Itel 2 + kl2| 2 + kl4|) 2 

O J — 00 

- - dx (\q 15 \ 2 + \q 25 \ 2 + ql 2 Qu + Q12Q14:) 2 ■ ( 68 ) 
a J -00 



6.3. Example of Z 6 -reduction. 

Let us impose the following Z 6 -reduction: 

B~\rf{x, t, X*))B = U{x, t, A), U(x, t, X) = Q(x, t) - XJ 

B = U7 ei _ e2 • U> e2 _ e3 ■ W e:j — e4 , -W ei — e5 , -W e5 — efj , B = M, 

where w ei - ej are the Weyl reflection with respect to the roots ej — ej of the so (12)- 
algebra. Then B(J) = J and the corresponding reduced potential Q red G so(12) 
takes the form: 

n rcd _ ( q \ 

Q ~U /' 
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where 



q = 
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Thus one derives the following 3 component MNLS system related to D$- algebra 
(here the independent fields are: (712(2;, t), qi3(x, t), and qu(x, t) : 

. dqu <9 2 9l2 „ ,2 . | ,2 , ol |2\ 

ia ~Of + ~q^2~ ~ 2 9i2(3|gi 2 r + 2|gi 3 | + 2 l9i4| ) 

+49i4(|9i3| 2 - I912I 2 ) + 2 9i 3 9i 2 - 2 9i 2 9* 4 ~ 2 9i 3 9i4 ~ 2 9i 4 9i 2 = 

. <9<?i3 d 2 qi 3 2 2 2 

ia ~Qf + -q^2~ ~ 2 9i3(3|gi 3 | +2|gi 2 | +2|gi 4 | ) 

+ 2 9i 2 9i 3 + 2 9i 4 9*3 + 4 9i39i29i4 + 4 9i39i49i2 - 4 9i29i49*3 = ( 69 ) 

. #914 . d 2 q U 2 2 , a I |2\ 

la ~gf + -Q^T ~ 2 9l4(|9l4| + 4 l9l2| + 4 l9l3| ) 

+ 4 9i2( 2 l9is| 2 - |912| 2 ) + 4 9i 2 3 9i4 - 4 9 2 3 9i2 - 4 9i 2 2 9i* 4 = 

The Hamiltonian of (69) is obtained from the general expression (33) by imposing 
the reduction constraints: 

H = flkin + -f^int, 

ffkin .ir* + 4=^ + 2 ^M) , (V0) 

a J -00 \ ox ox ox ox ox ox J 

00 



tfint = - r dx (6|gi 2 | 4 + 6|gi 3 | 4 + |9i4| 4 + 8|gi 4 | 2 (|9i2| 2 + |9i3| 2 )) 

O J—oo ^ ' 

6 f°° 

+ - dx (q 2 2 qtt + q\ A q\l - q 2 3 q{t - g? 4 g*§ - g 2 2 9i! - 9 2 3 9ii), 

O J—oo 

12 f°° 

H / dx ((|<?i 2 | 2 - 2|gi 3 | 2 )(gi 2 q* lA + q u g* 2 )), 

O J—oo 



12 



H / dx (9i29i49i3 + 9i39i49i 2 )), 

O J—oo 



The soliton solutions of the reduced MNLS require additional efforts. The prob- 
lem is that the generic expression for the dressing factor (57) does not satisfy the 
reduction conditions. 
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7. Reductions and the scattering data 

The reduction conditions (63) and (64) imposed on the potential of the Lax operator 
(2) induce an invariance conditions for the corresponding fundamental analytical 
solutions 

Typel: B-\ x + (x, t, X*)?B = ( C ^ A) a _° (A) ) ( x ~(x, t, A))" 1 , (71) 

Type II: C~ 1 (x ± (x, t, X*))*C = x ± (x, t, A), (72) 
and for the scattering matrix (24) 

Typel: B^T^t, X*)B = {T{t, X))' 1 , (73) 
Type II: C _1 T*(i, A*)C = T(t, A). (74) 
If we represent the internal automorphisms S a of type I in the form 

_(B + \ 

the reduction conditions on the matrix blocks a ± (A), b ± (X), cr^A) and d ± (X) in 
the scattering matrix and its inverse (24) read: 

B + (a + (X*))^B^ = c+(A), S + (6+(A*))ti? I 1 = d~(X), (75) 

B_(6-(A*))tfi+ 1 = d+(A), J B_(a-(A*)) t 5l 1 = c"(A). 

For the other reductions of type II the internal automorphisms preserving J up to 
a sign has off-diagonal block structure as follows: 

s - ( C A 

and the matrix blocks of T(A) are constrained by 

C+(a-(A*))*C+ 1 = a+(A), C_(a + (A*))*Cl 1 = a~(A), (76) 
C+^+CA*))*^ 1 = -6"(A), C_(6"(A*))*C+ 1 = -6 + (A). 

Similar reduction constraints could be written also for the inverse T(A) to the scat- 
tering matrix T(A). 

Our last remark here is that the reductions described above can not be applied to 
any generic NLEE. Indeed, eq. (29) will be compatible with the reduction only if 
the dispersion law /(A) satisfies 

BZ 1 (/+(A*))t B + = fl\X) BZ 1 (/-(A*)) t BZ 1 = fz\X). (77) 

Any generic NLEE is compatible with reduction of Type II if the dispersion law 
complies with 

C + (/-(A*))* CZ 1 = f + (X) C- (/ + (A*))* CZ 1 = /-(A). (78) 
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8. Conclusions 

We have described new systems of MNLS type obtained as Z 4 and Zg-reductions 
of the MNLS related to a D.III type symmetric space. The Hamiltonian formalism 
and the theory of A-operators for MNLS related to the relevant simple Lie algebras 
are briefly discussed. We show how the method, presented in [23] for the N- 
wave equations and their gauge equivalent systems can be extended to MNLS type 
systems [25]. The reduction of the multi-component nonlinear Schrodinger (NLS) 
equations on symmetric space C.I ~ Sp(2p)/U(p) for p = 2 is related to spinor 
model of Bose-Einstein condensate. Other interesting reductions of MNLS type 
equations were reported in [9] and a systematic study of the problem is on the way. 
These results can be extended and the reductions of MNLS-type equations related 
to other symmetric and homogeneous spaces can be explored. As a result one 
can systematically obtain and classify new integrable systems of MNLS type. The 
method is explicitly gauge covariant and can also be applied to their gauge equiv- 
alent systems of Heisenberg ferromagnet type. Such research would entail a volu- 
minous calculations and will be continued in subsequent publications. 
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